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ABSTRACT 

Fullerenes have recently been identified in space and they may play a significant role in the gas and 
dust budget of various astrophysical objects including planetary nebulae (PNe), reflection nebulae 
(RNe) and Hn regions. The tenuous nature of the gas in these environments precludes the formation of 
fullerene materials following known vaporization or combustion synthesis routes even on astronomical 
timescales. We have studied the processing of hydrogenated amorphous carbon (a-C:H or HAC) nano- 
particles and their specific derivative structures, which we name "arophatics" , in the circumstellar 
environments of young, carbon-rich PNe. We find that UV- irradiation of such particles can result in 
the formation of fullerenes, consistent with the known physical conditions in PNe and with available 
timescales. 

Subject headings: circumstellar matter - infrared: general - ISM: molecules - planetary nebulae: 
general - stars: AGB and post-AGB 



1. INTRODUCTION 

The gas and dust that is ejected by dying stars and 
from which new stars will form, is constantl y being 
proce ssed chemically as well as physically (e.g., Tielens 
2005). A particularly rich and diverse route is followed 
by carbon. In the cool surroundings of carbon stars 
- evolved stars that have dredged up large amounts 
of carbon to the stellar surface - more than 60 indi- 
vidual molecular species have been id entified, including 
benzene, polyynes an d 
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cyanopolyy nes ( Kwok 20091 |Cer- 
Pardo et all 120070. In addition, 
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the spectra reveal diverse dusty materials such as hy- 
drogenated amorphous carbon (HAC or a-C:H) and sil- 
icon carbide (SiC). As the mass- loss process strips the 
outer layers of the star, a hot (T > 30,000 K) central 
white dwarf becomes exposed, whose strong UV irradia- 
tion further processes the ejecta and makes them visible 
as a planetary nebula (PN). A key component of the car- 
bon dust inventory in PNe are polycyclic aromatic hy- 
drocarbons (PAHs) and PAH-like species, a class of large 
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and hardy carbon aceous species whose for mation mecha- 
nisms are unclear ( Cherchneff et al.|2000 ). Although not 
a single PAH member has been identified in space, their 
characteristic spectral features at infrared wavelengths 
are observed throughout the Universe, from which it is 
inferred that they reside ubiquitously in space and c ould 
lock up as much as 15% of the cosmic carbon (e.g., Tie- 
lens|[2005 ) 



Two other large aromatic species have been recently 
identified in space: Ceo and C70. These are the best- 
known members of the family of fullerenes, a class of 
molecules made of hexagonal and pentagonal aromatic 
carbon rings, fused in the shape of a hollow sphere or el- 
lipsoid. The most abundant of these molecules, the Buck- 
minsterfullerene Ceo, has the structure of an old- fashion 
black and white soccer ball. The near-spherical carbon 
configuration of these two molecules, which have a closed 
surface, a closed-shell electron distribution and an al- 
most unstrained network, results in a very high stabil- 
ity against dissociation and prolonged exposure to high 
temperature. Fullerenes, in particular Ceo, have peculiar 
and appealing photochemical, electrochemical and phys- 
ical properties, which can be exploited in various fields, 
from nanotechnology to medicine, to space science. 
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Ceo and C70 were discovered during laser ablation ex- 
periments on graphite targets, ai ming to study long car- 
bon chains in interstellar clouds ( |Kroto et al.|fl 985). Be- 
cause of their remarkable stability, fullerenes appeared 
particularly suited to survive the harsh conditions of the 
interstellar medium. Their unique properties indicate 
that they may play an important role in organic astro- 
chemistry and astrobiology. Cosmic fullerenes remained 
elusive until the recent di scovery of C^p and C70 in the 
planetary nebula Tc 1 by|Cami et al. (2010). Fullerenes 



have since b een confirmed m many more evolved star en- 
vironments ([G arcia-Hernande z et al.|2010[|2011b[|Z hang 
Kwok||20lil JCIayton et al.||Mll IQielen et alH MH 



Roberts et al.N2012| JEvans et al.|J2012| , as well as m 
young stella r objects ([ Roberts et al. 20 l2p, reflection neb- 
ulae (RNe) flSellgren eTal.12010] IP eeters et al.|2012| ) and 
photodissociation regions ( |Rubin et al.|2UllJ ). 

Fullerenes can be efficiently synthesized in the labo- 
rato ry by the vaporization of ca rbon rods in an electric 
arc (Kratschmer et al. 1990a|b[ ) and fro m hydrocarbon 
combustion under o ptimised conditions (Howard et al. 
1991 Nano-C||2004[ ). However, the formation routes of 
fullerenes m space are still unknown. We review the cur- 
rently known formation mechanisms of fullerenes in the 
laboratory (§|2]), showing why these methods would not 
work in space (§ [3). We propose an alternative path- 
way, consistent with astrophysical conditions, based on 
the photo-processing of a family of carbonaceous species 
which we name "arophatics" (§]4|. We discuss the ob- 
servations of fullerenes in PNe (spf and finally we sum- 
marize our conclusions (§[6]). 

2. FULLERENE FORMATION ON EARTH 

2.1. The discovery of fullerene 

The original discovery of C bq and C70 dates ba ck to 
1985 with the experiment of Kroto et al.| ([1985). In 
this experiment, carbon species were vaporized from 
a graphite target into a He flow with tunable pres- 
sure/density. The vaporization was done by a pulsed 
laser, the resulting carbon clusters were expanded in a su- 
personic molecular beam, photoionized and their masses 
measured by a time-of-flight mass spectrometer. Vapor- 
ization at low He pressure (less than 10 torr) leads to 
a broad distribution of clusters, with 38 to 120 carbon 
atoms (always even numbers), with Ceo and C70 present 
but not the dominant species. At a pressure of 760 torr 
the C60 and C70 peaks clearly dominated. When the ex- 
pansion of the He+carbon cluster mixture was delayed, 
the resulting mass distribution was completely domi- 
nated by Ceo - 50% of total large cluster abundance, and 
C70, with 5 % of total large cluster abundance. These 
distributions have been interpreted as due to the increas- 
ing number of He-cluster and cluster-cluster collisions, 
resulting in the "survival of the fittest", i.e. the more 
stable clusters. The famous icosahedral soccer ball struc- 
ture was proposed for the first time for Ceo , and later 



confirmed by spectroscopic studies ( Kratschmer et al i 
1990a|b Taylor et al. 1990). However, the elementary 
reaction mechanisms occurring during the "thermaliza- 
tion" process remained unknown. 

2.2. The role of temperature: mutually exclusive 
formation of PAHs and fullerenes 



Jager et al. (2009) experimentally studied the gas- 
phase formation of carbonaceous compounds. Their first 
experiment was essentially a replication of the Kroto ex- 
periment, i.e. laser vaporization of a graphite target in 
a quenching atmosphere of He or He/H2 with pressures 
between 3.3 and 26.7 mbar. The vibrational temperature 
of the laser-induced plasma was estimated to be between 
4000 and 6000 K for the laser power densities used. In 
a second set of experiments, the laser-induced pyroly- 
sis of C2H2 and CeH6 using a pulsed laser with high 
power densities was studied. The resulting condensa- 
tion temperatures were above 3500 K, comparable with 
the vaporization experiment. In a third set of experi- 
ments, laser-induced pyrolysis of the same hydrocarbon 
precursors was carried out, but this time at a much lower 
temperature (max 1700 K). 

The analysis of the condensation products showed a 
striking effect of the temperature. During the high 
temperature experiments, only fullerene-like soot and 
fullerenes were produced. During the low temperature 
experiments, only soot and PAHs were formed (100 % 
PAHs at ~ 1000 K). The absence of PAHs in the graphite 
vaporization experiment is maybe not so surprising, due 
to the lack of hydrogen, but the results of the pyrolysis 
experiments clearly tell that the temperature determines 
which kind of condensates will be formed. Moreover, the 
two pathways seem to be mutually exclusive, at least un- 
der these experimental conditions: fullerenes and PAHs 
cannot be formed together. It is important to note that 
high pressure has been used during the experiments to 
concentrate the condensation within a small volume, i.e. 
to mantain a high density of the condensing species. 

2.3. Molecular dynamics simulations of fullerene 
formation 

The laboratory experiments described in the previ- 
ous sections do not tell us how fullerene is formed from 
graphite and hydrocarbons, but only under which condi- 
tions this happens. The elementary reaction mechanisms 
involved in the formation of fullerenes can be investi- 
gated perfor ming quantum chemical molecular dynamics 
simulations ( | Zheng et al.|[2005| |Irle et aT]|2006 Zheng 
et al.||2007|). In one set 01 simulations a hot carbon va- 



por was reproduced by putting 40 ra ndomly oriented C 2 
molecules in a tiny 20 A cubic box (Zheng et al. |2005 ). 
The system was let to evolve, and extra C2 molecules 
were added at fixed times in order to simulate an open 
environment. During the first stage, the system was kept 
at a constant temperature of 2000 K and giant carbon 
cages self-assembled. During the second stage of the sim- 
ulation the temperature was raised to 3000 K, and the 
shrinking of these hot giant fullerenes down to C65 was 
observed. 

The formation of the giant cages starts with nucleation 
of polycyclic structures (hexagons and pentagons) from 
entangled polyyne chains. This is followed by growth of 
the structure and finally cage closure similar to the self- 
capping of open-ended single- walled nanotubes. Because 
of the rapid gain of energy due to cage closure, the giant 
fullerenic cages are produced in a vibrationally highly 
excited state. The excess energy has to be dissipated, ei- 
ther by unimolecular dissociation or collision with other 
carbon clusters or carrier gas atoms. The simulations 
show that the newly produced giant fullerenes inevitably 
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shrink to smaller sizes. All the road maps to fullerene 
formation proposed by previous models were associated 
with intermediate structures that are in thermodynamic 
equilibrium, while a hot carbon vapor is indeed a system 
far from thermodynamic equilibrium. In the model de- 
veloped by Irle, Morokuma and collaborators, for the first 
time the dynamic self-assembly of fullerene molecules oc- 
curs as an irreversible process emerging naturally under 
the nonequilibrium conditions typical of a hot carbon va- 
por. 

2.4. Direct formation of fullerenes from graphene 



Chuvilin et al. (2010) have shown that fullerenes can 
form directly from graphene. In their experiment, a sheet 
of graphene was exposed to a beam of energetic 80 keV 
electrons. The energy transferred to the graphene re- 
sults in the loss of carbon atoms from the edges of the 
sheet. If the size of the graphene flake is not too big 
(less than a few hundreds of carbon atoms), the loss of 
carbon atoms at the edge will lead to the formation of 
pentagons, which triggers the curving of graphene into 
a bowl-shaped structure. Further carbon removal from 
the edges using the electron beam will reduce the size of 
the curved structure until it is sufficiently small to zip 
up its open edges and isomerize into a closed fullerene 
structure. 

2.5. Formation of fullerenes via closed network growth 



Dunk et al.| (120121) claim having identified the funda- 
mental processes leading to the formation of fullerenes in 
a recent experimental work, based on laser vaporization 
of pure Ceo in the presence of carbon vapor. Accord- 
ing to this study, fullerenes would self-assemble bottom- 
up through a closed network growth (CNG) mechanism 
based on the ingestion of atomic carbon a nd C2 clus- 
ters. It should be noted that the work from lDunk et al.l 
( |2012| ) provides experimental evidence of the growth of 
larger fullerenes from pre-existing Ceo only. Because the 
experiments show that the CNG of larger fullerenes does 
not result in the production of Ceo, it is deduced that 
Ceo formation must be a result of CNG from smaller 
fullerenes. However, the initial formation mechanism of 
such smaller fullerenes is still under debate. 

3. FULLERENE FORMATION IN SPACE: WHY 
TERRESTRIAL METHODS DO NOT WORK 

Although all the chemical ingredients and the required 
temperatures for graphite vaporization and hydrocarbon 
combustion can be found in astrophysical environments, 
the densities are far too low to proceed to fullerene for- 
mation on reasonable timescales. In space, the most 
favourable conditions for fullerene formation via vapor- 
ization/combustion are found in a post-shock gas. In 
such environments the fullerene building blocks, i.e. C2 
groups and poly cyclic species coming fro m the fragmen- 
tation of PAHs ( jMicelotta et al.j2010a|b ) , ought to exist 
in the required vibrationally excited states. However, the 
post-shock carbon densities are low and represent a seri- 
ous obstacle to fullerene formation. We derive a scaling 
rule relating the timescale f for fullerene condensation 
with the initial density n of the carbon gas. This scaling 
rule is then used to estimate the time required for the co- 
agulation of fullerenes from a carbon gas with interstellar 
densities (n ~ 10 -4 carbon cm -3 ). 



Consider an ensemble of N p particles that can occupy 
N v cells. The probability of having one particle in one 
cell is p = N p /N v <C 1. The probability P of having two 
particles in the same cell at the same time is described 
by the Poisson statistics: 



P(2;p) 



P 



(1) 



We have that N v = V/6v, where V is the total volume 
and Sv is the volume of a single cell. We can then rewrite 
the single occupancy probability as p = (N p /V) 5v = 
n Sv. 

The probability P of having interaction between two 
particles at time N t can be written as: 



P(int;JVt) = (l-A) 



Nt-l 



A 



(2) 



where N t — t/r, t being the time and r the characteristic 
time of a single interaction (depending on the nature of 
the interaction itself). The term A represents the prob- 
ability of interaction between two particles at a generic 
time t. The probability of having interaction between 
two particles up to the time N t is then given by the fol- 
lowing equation: 

Nt-l 

P(int; < N t ) = ^ (1 - A) fe A = 1 - (1 - X) Nt (3) 

k=0 

If A is small, we can rewrite Eq. [3] as: 

P(int; < N t ) ~ 1 - e~ XNt ~ 1 - e~ x ( * /r) (4) 

The probability A can be written in the following way: 

\ = AP(2-p) ~A(nSv) 2 (5) 

The term A includes the characteristics intrinsic to the 
interaction (which we do not need to specify) while 
P(2;p) is the probability from Eq.fl] The approximation 
resulting in the right-hand term 01 the equation is valid 
for p — )> 0. Combining Eqs. [4]and[5j with ASv 2 = B, we 
obtain 

P(int; < Nt) ~ 1 - e A ^ 2 ^ ~ 1 - e BnH ' T (6) 

From Eq. [6] we can derive the timescale for condensa- 
tion: 

Because both r and B are independent from the initial 
particle density n, we can write the following scaling rule 



T_2_ 



r 2 



(8) 



where f\ and f 2 are the condensation timescales corre- 
sponding to densities n\ and n 2 respectively. 

Knowing the density and timescale in simulations, n\ 
and fi, and the density of the interstellar gas, 77,2, we can 
then calculate the fullerene condensation time in space, 
f 2 . We adopt n\ = 10 22 carbon cm -3 and f\ = 50 ps 
( |Zheng et al.||2005[ ). The hydrogen density of the post- 
shock gas is ~ 1 cm -3 , but carbon is 10 -4 less abun- 
dant, so n 2 = 10 -4 carbon cm -3 . From Eq. |8] we ob- 
tain T2 = 1.6 x 10 34 yr, which is longer than the age of 
the Universe. We also estimate the carbon density n 2 
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required to form fullerenes in a shocked gas before the 
temperature drops below 1000 K. At this temperature 
the condensation of fullerenes cannot occur. For a typical 
interstellar shock with a velocity of 125 km s _1 the cool- 
ing time (after which the temperature of the post-shock 
gas drops below 1000 K) is about 2xl0 5 yr. Under these 
conditions, the required carbon atom density for the con- 
densation of fullerenes is = 3xl0 10 carbon cm -3 , and 
in space such densities do not exist outside of stellar at- 
mospheres^ 



Berne & Tielens (2012) propose a scheme for the di- 
rect formation of (Jqo from PAHs in the reflection nebula 
NGC 7023 based on grap hene dissociation experiments 
fromlChuvilin et al.|poTo|. In NGC 7023 fullerenes and 



PAHs appear to coexist m the same location close to the 
central star, while only PAHs are detected further away 
from the star. The absorption of UV photons from the 
central star fi rst completely dehydrogenates the PAHs 
(Tielens 2005), producing graphene flakes. Further pho- 
ton absorption induces the loss of carbons f rom the edges 
of the graphene sheet. This, according to |Chuvilin et al.] 
( |2010| ), gives rise to pentagonal defects which will induce 
the curvature of the sheet and the subsequent closure of 
the cage to form a Ceo molecule. 

There are some issues related to the applicability of the 
results from such experiments to interstellar conditions. 
First of all, it is not possible in space to tune the number 
of carbons removed from the graphene flake in order to 
get the right size required for cage closure. The number 
of ejected carbons depends on the photodissociation pro- 
cess, which is rather slow. To form fullerenes during the 
lifetime of NGC 7023, the formation mechanism needs 
to be fine-tuned. The precursor PAH molecule needs to 
contain about 70 carbon atoms in order to be able to 
form Ceo, and the most favourable dissociation condi- 
tions are required. The carbon atom loss is governed by 
the parameter Eq, which describes the dissociation rate 
of a highly exc ited PAH molecule using an Arrhenius law 
Tielens|2 005). The value of Eq for PAH dissociation un- 
er interstellar conditions is not well established. For this 



reason, a ran ge of values has to be considered (Micelotta 



et al.||2010b). In the scheme proposed by Berne fc He- 



lens (12012), a low value has to be chosen, ~EJq = 3.65 



eV. IncT eed, if the dissociation parameter is higher, e.g., 
4.6 eV ( jMicelotta et al.||2010b[ ), then the fragmentation 
process is far too slow and the proposed conclusion no 
longer holds: PAHs with 70 carbon atoms require longer 
than the lifetime of NGC 7023 to fragment, therefore 
they cannot be the precursors of the Ceo observed at the 
present epoch. 

Thus, none of these mechanisms are feasible in the as- 
trophysical settings where fullerenes have been detected, 
and hence, alternative routes are needed. 

4. OUR PROPOSED FORMATION MECHANISM: 
FROM A-C:H TO FULLERENE THROUGH 
AROPHATIC CLUSTERS 

Trace amounts of fullerenes have been identified among 
the dissociation products of lase r- irradiated HACs (|Scott| 



et al.|1997||Gadallah et al.|2011| ). Together with fullerene 
emission, the three PNe Tc 1, SMP LMC 56 and SMP 
SMC 16 show broad emission plateaus between 6-9 /im 
and between 10 - 13 /im, compatible with emission from 
H-rich amorphous carbon nanoparticles ( |Bernard-Salas| 



et al. [2012 ). This suggests that a viable formation route 
m sp ace could exist that starts from HA C processing 
(e.ff., |Garcfa-Hernandez et al.||2010[|2011aD . 

In this study we locus on the formation of fullerenes in 
the circumstellar environments of proto-planetary nebu- 
lae (PPNe) and PNe. HACs are an abundant carbona- 



ceous dust component in both these types of objects (Gr 
ishko et al ||2001| |Goto et al l|2003| ). In addition, oneot 



the main characteristics of PNe is their strong radiation 
field, which will process the HACs. 

4.1. Arophatic clusters 



In t heir experiments, |Smit"h (|1984 ) and Scott & Duley 
(1996) showed that HA(J/a-C:H decomposition leads to 



a low density material (p < 1.5 g cm 3 ) . In p articular, 
|Scott fc Duleyl fll996| ) and [Scott et al.] fll997| ) showed 
that the decomposition end-product is an aerogel-like 
material consisting of weakly-connected aromatic "proto- 
graphitic" or PAH-like clusters in a "friable network". 
The dehydrogenation of small carbonaceous particles, 
containing less than a few hundred C atoms, will likely 
lead to structures that resemble the particles seen by 
|Scott et al.| (1997) in their experiments and/or the lo- 
cally aro matic polycyclic h ydrocarbons (LAPHs) pro- 



posed by Pe trie et al. (2003]). The particles in the [Scott 
et al. ( 1997p experiments are found to be small alka- 



nes, unsaturated carbon-chain radicals and small dehy- 
drogenated PAH- like species. After the release of these 

1997[) observed a strong 



lighter molecules |Scott et al. 

mass peak near 500 amu [~ 40 C atom, proto-graphitic, 
aromatic clusters) having IR spectra similar to that of 
HAC/a-C(:H) with absorption and emission features at 
3.3, 3.4 and 6.2 /im. The results o f these experiments 
and recent modelling ( |Jones|[2012c[ ) therefore imply an 
intimate mix of aromatic clusters and aliphatic groups 
in small carbonaceous species. We note th at the LAPH 
structures studied by |Petrie et al.l (|2003[ ) have compo- 
sitions in the range C19H22 to C36H32 and contain aro- 
matic domains with only a few rings per domain, which 
are linked by bridging aliphatic ring systems, and should 
have IR spec tral properties s imilar to HAC/a-C(:H) ma- 
terials (e.g., Jones|[2012a|b|c ) . 

The conditions required lor HAC dehydrogenation are 
common in space, we can thus expect the formation of 
LAPH-like species where the aromatic components are 
bridged by olefinic carbon-containing rings and chains. 
The key component in these structures is the presence of 
the aliphatic bridging groups linking the aromatic clus- 
ters, in what we here call "arophatic" structures. Such 
structures, analo gous to an aliphatic per iphery on PAH 
molecules (e.g., |Rauls fc Hornekger|[2008| ) , could p rovide 
sites fo r molecular hydrogen formation m PDRs ( |Jones| 
|2012b[ ) . Size-dependent tt — 7r* band considerations indi- 
cate that the intrinsic aromatic clusters are indeed most 
likely to consist of 'isolated' two- and three-ring sy stems 
in such small hydrocarbon particles ( Jones||2012c| ). The 
non-aromatic carbon atoms are in the form of "more 
flexible" short olefinic and/or aliphatic bridging struc- 
tures such as -CH=CH- and -CH2-C H2-. This generic 
type o f structure was also proposed by Gr ishko fc Duley] 
( 2000 ) , who suggested that aliphatic cyclo-hydrocarbons 
are the precursors to aromatic cluster formation in HAC 
materials. The aliphatic bridging structures must be 
quasi-orthogonal to the aromatic species that they bridge 
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and will be unable to transform into aromatic structures 
without compromising (i.e., destroying) the entire struc- 
ture. These active sites could, likely, interchange between 
sp 2 olefinic and radical sp 3 states with dangling bonds, 
in a kind of isomerisation "flipping" process, which could 
play a key role in cage closure and fullerene formation in 
the solid state; 

The work of Ferrari fc Robertson| (2QQ4) indicates that 
small clusters within a-C(:H) materials are chain-like, 
while the larger cluste rs are three-dimension al, cage-like, 
sp 2 structures and as Chuvilin et al.| ( |2010| ) show these 
cage-like structures have a natural tendency to close-up 
and form fullerene-like and fullerene structures. 

The "arophatic" clusters are structures evolving from 
larger, amorphous hydrocarbon particles and, from a 
topological point of view, they will inevitably end up 
with a tube-like, cup-like or cage-like appe arance (e.g., 
see th e two right hand species in Fig. 1 of |Petrie et al.| 
2003). Such structures are 'curled- up' or folded-over', 



delected graphite (DG) or graphene sheets ( |Jones 2012a), 
and we can therefore represent t hem by their ' tiattened- 
out' equivalents as per Fig. 12 in | Jones (2012a) but with 
aromatic domains reduced in size to only a lew rings per 
aromatic cluster. An example of such an un-wrapped 
"arophatic" is shown in Fig. [I] while Figs. [2] and [3] show 
the folded-around version of this size-limited structure. 

The aromatisation of the aliphatic bridges to form 
pentagons can provide a viable route to fullerene forma- 
tion in the solid phase via these "arophatic" species. The 
key process is the dehydrogenation of the aliphatic and 
olefinic "bridges" which naturally leads to aromatic pen- 
tagon formation, as shown in Fig. H] for small a-C:H par- 
ticle fragments. It appears that there is a link between 
small carbonaceous particle dehydrogenation in intense 
UV radiation fields and the formation of fullerenes, which 
could explain the observation of fullerenes in circum- 
stella r regions associated with small hydrocarbon grains 



e.g., 



Cami et al. 2010 



2011b 



Cami et al 2011 



Garcia-Hernandez et al. 2010 



Bernard- S alas et al. 2 012[ )T 



4.2. Fullerene formation in space 

The mechanism that we propose for fullerene for- 
mation in space is the following. The starting point 
is large, nm-sized, H-rich a-C:H particles containing 
more than a hundred carbon atoms. UV photolysis 
by photons with energy above 10 eV dehydrogenates 
these grains ( Jones|[2012b[ ) and leads to further changes 




Figure 1. An example of an un-wrapped small "arophatic" cluster 
(C46H38). The filled circles represent carbon and the dotted bonds 
indicate the structure-wrapping connections. The folded-around 
version of this structure is shown in Fig. [2] 




Figure 2. The wrapped version of the small "arophatic" cluster 
(C46H38) shown in FigJT] Carbon is represented in dark grey, 
hydrogen in light gray. The color version of this figure (carbon in 
red, hydrogen in gray) is available in the electronic edition of the 
Astrophysical Journal. 




Figure 3. A snap-shot from Movie 1, which illustrates the 3-D 
structure of the small "arophatic" cluster C46H38 schematically 
represented in Fig. [2] Carbon is in dark grey, hydrogen in light 
gray. Movie 1, with carbon in red and hydrogen in gray, is available 
in the electronic edition of the Astrophysical Journal. 



in the molecular structure - a process that we call 
STIR-processing (Structural Transformations in Intense 
Radiation fields). STIR processing requires constant 
heating (by photon absorption) to keep the particle in 
a vibrationally excited state that promotes structural 
transformations leading to the closing of cage-like struc- 
tures within the a-C:H nanoparticle. Amorphous carbon 
nanoparticles can be maintained at thermal-equilibrium 
temperatures of the order of 100-150 K by photon ab- 
sorption (for H-rich a-C:H particles with band gaps larger 
than 1 eV) even at thousands of astronomical units from 
the central stars of PNe, due to their low emissivity 
at the long waveleng ths where they radiate thermally 
(Bernard-Salas et al.|2012). However, at these same dis- 



tances, almost completely dehydrogenated nanoparticles 
have temperatures of only tens of Kelvin, as expected. 

In order to form fullerenes, the arophatic clusters re- 
sulting from STIR processing must be larger than the 
particle previously shown (C46H38), and contain more 
than a hundred carbon atoms, as exemplified by the clus- 
ter in Figs. [5] and [6] 
A key step in the route to fullerenes is the formation 
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bon is represented in dark grey, hydrogen in light gray. The color 
version of this figure (carbon in red, hydrogen in gray) is available 
in the electronic edition of the Astrophysical Journal. 




Figure 4. A schematic view of three possible routes to aromatic 
pentagon formation (labeled A, B and C), leading to a loss of pla- 
narity, in nm-sized a-C:H particles via aliphatic and olefmic ring 
system dehydrogenation. The filled circles represent carbon. 

of 5-membered rings. This can occur via dehydrogena- 
tion during the STIR processing of the parent a-C:H or 
of the emerging arophatic cluster, in a similar fashion 
as shown in Fig. [4] for small a-C:H particles. Dehydro- 
genation is a viable way to introd uce 5-membered rings , 
as indicated by the experiment of |Gadallah et al" (2011) 
where hydrogenated amorphous carbon material is UV 
irradiated simulating the interaction of cosmic carbona- 
ceous nanoparticles with UV radiation in space. They 
typically find strongly distorted graphene layer struc- 
tures, small fullerene-like carbon structures and aro- 
matic clusters containing 3 — 8 rings (with diameters 
c± 0.4 — 0.7 nm). This clearly shows that UV irradia- 
tion not only causes dehydrogenation but also forces the 
structure to curl up because of the introduction of pen- 
tagonal rings. In fact, the UV-induced dehydrogenation 
is the cause of the formation of 5-membered rings. Thus, 
in a single step, dehydrogenation fulfills two essential 
conditions for the formation of fullerene molecules: 1) 
removing hydrogen atoms and 2) triggering pentagonal 
ring formation in the evolving structure. The fullerene 
precursor is a big and complex arophatic cluster, with the 
aromatic domains being non-co-planar. The 3-D struc- 




Figure 6. A snap-shot from Movie 2, which illustrates the 3-D 
structure of the large and complex "arophatic" cluster C104H80, 
precursor of cosmic fullerenes and schematically represented in 
Fig. [5] Carbon is in dark grey, hydrogen in light gray. Movie 2, 
with carbon in red and hydrogen in gray, is available in the elec- 
tronic edition of the Astrophysical Journal. 



ture favors the formation of 5-membered aromatic rings, 
and then curvature, even in such a big structure, via 
photo-induced dehydrogenation. The particle emerging 
from the dehydrogenation process is a giant fullerenic 
cage similar to the ones formed in simulations ( |Irle et ah] 
|2006| . The initial, STIR-driven (catastrophic) photolysis 
and the following cage closure induced by the formation 
of pentagonal rings will result in a highly excited struc- 
ture. The cage will release the excess energy via ejection 
of C2 molecules, shrinking down to the islands of stabil- 
ity represented by Ceo and C70 ( |Irle et aL]|2006| ). The 
down-sizing occurs through unimolecular decomposition 
of the vibrationally excited cage, resulting in the ejec- 
tion of carbon molecules. To model the unimolecular 



dissociation we adopt an Arrhenius law (Tielens 2005) 



which describes the evolution of the dissociation rate as 
a function of the temperature of the dissociating species: 



k s (T) = k exp(-^oAT) 



(9) 



Cosmic Fullerenes from Arophatic Clusters 
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Table 1 

Timescale, t s , for shrinking of the fullerenic cage C144 to C60- 



T 


k s 


t s 


(K) 


(C/s) 


(s) 


1000 


1.8xl0 9 


4.7xl0- 8 


300 


1.1x10 s 


0.7xl0- 3 


200 


l.OxlO 2 


8.4X10- 1 


100 


l.OxlO- 7 


8.4xl0 8 



Note. — T is the temperature of the cage and k s is the rate for 
carbon ejection measured in carbon atom lost per second. 

The quantity fc s is the rate for carbon ejection measured 
in carbon atom lost per second (C/s), T the vibrational 
temperature of the fullerenic cage and fc the Boltzmann 
constant. The unimolecular decomposition is governed 
by two parameters: the pre-exponential factor fco and the 
dissociation energy Eq. To estimate the value of these 
two parameters, we use two sets of simulations which 
follow the shrinking of giant fullerenic cages at two dif- 
ferent temperatures, 2000 and 3000 K, and provide the 
corresponding dissociation rates. The average dissocia- 
tion rate k s is 1.5 xlO 10 C/s and 3 xlQ 1Q C/s for T = 2000 
K and T = 3000 K respectively fl Zheng et aLl|2007| ). We 
obtain the following values for fco and ~E$: 

fc = 1.2 x 10 11 C/s 
£ = 0.36 eV 

The time t s required for a big fullerenic cage to shrink 
down to Ceo is given by: 



t B (T) = N ej /k s (T) 



(10) 



where 7V e j is the number of carbons that have to be 
ejected from the cage and fc s is the ejection rate from 
Eq. [9] We consider a fullerenic cage containing 144 car- 
bon atoms, i.e. with size comparab le to the size of th e 
particles considered in simulations (Zheng et al. 2007), 
therefore, 7V e j = 84. The shrinking timescales tor vari- 
ous temperatures of the initial cage are reported in Ta- 
ble [T] This process requires relatively short time scales: 
to shrink a large (144 carbon atoms) fullerene cage to 
Ceo requires only about 0.8 s at 200 K, and increases to 
27 years at 100 K. 

Ceo and C70 are the only, and therefore smallest, C n 
fullerenes with n < 70 satisfying the Isolated Pentagon 
Rule (IPR). This implies (almost) spherical, thermo- 
dynamically favored, closed shell structures particularly 
suited to survive harsh conditions. The shrinking stops 
at Ceo because of the very high Arrhenius dissociation en- 
ergy Eq for C2 elimination to C58, Eq ~ 10 eV. For com- 
parison, the dissoci ation energy for Cb 2 decomposition 
to C 60 is E ~ 7 eV flTomita et al. 12001 j ). The estimated 
temperatures of a few hundre ds 01 Kelvin for fullerenes 



2010 



Bernard- 



detected in plane tary nebulae ( Cami et al i 
Salas et al.|2012 ) are far from the approximately 1 UU0 K 
required to overcome this decomposition barrier ( |Leifer 
et"aI1[T995| . 

Our model proposes the formation of a single Ceo (C70) 
molecule from a single large arophatic cluster. This im- 
plies that the parent a-C:H particles must be nm-sized, 
and therefore contain a hundred or more carbon atoms, 
because of the dissociative loss of carbon (and hydro- 
gen) atoms that will occur during structural transforma- 



tion. Note that the large arophatic clusters can be fur- 
ther fragmented by the same radiation field that caused 
their formation from the parent a-C:H grain. Most of 
the formed arophatics will fragment into structures other 
than fullerenes, resulting in a fullerene formation process 
that is rather inefficient, consistent with the low fullerene 
abundance observed in space so far. 

5. FULLERENE OBSERVATIONS IN PLANETARY 
NEBULAE 

The peculiarities of the PNe where fullerenes have 
been detected are their intense UV radiation fields, suffi- 
cient to sustain photo-chemical processing of HACs, and 
the ablation of carbonaceous material from the dense 
circumstellar matter formed during their AGB phase. 
Thus, in these objects abundant, nm-sized, H-rich, a- 
C:H dust is being liberated and undergoing catastrophic 
UV-photolysis, which leads to STIR processing. In Tc 1, 
Ceo and C70 have been detected far fo r the central star 
(8000 AU - iBernard-Salas et al l|2012| ), but even at this 
distance the U V irradiation is sufficient to keep the parti- 
cles in a vibrationally excited state (T ~ 100-130 K) and 
to trigger and sustain the STIR processing of the a-C:H 
dust and the shrinking of the resulting giant fullerenic 
cages down to Ceo and C70. 

It appears that the young PNe where fullerenes have 
been observed are the optimised factories for their forma- 
tion because their stellar emission peaks in the required 
10 eV photon energy regime. In their pre-PNe stages 
these objects do not have sufficient UV photons to dehy- 
drogenate and heat the dust, and in their later stages the 
radiation field is so much harder that dissociative losses 
will destroy the H-rich carbonaceous particles before they 
can re-structure to form fullerenes. 

6. CONCLUSIONS 

We studied the formation of fullerenes Ceo and C70 
in planetary nebulae, for which a satisfactory formation 
path is still missing. 

In space, the low density of the gas precludes the for- 
mation of fullerene materials following known vaporiza- 
tion or combustion synthesis routes even on astronomical 
timescales. 

The direct formation of Ceo through dissociation- 
induced curvature of dehydrogenated PAHs requires a 
very specific tuning of the dissociation parameters, and 
thus appears unlikely. 

The scheme for fullerene formation in space that we 
propose is based on a top-down approach. The top of 
the tree is represented by large a-C:H/HAC particles, 
i.e. materials whose presence in space has been firmly 
established. 

UV photolysis promotes structural transformations in 
the parent a-C:H particle with the formation of what we 
call "arophatic" clusters: 3-D, hollow structures char- 
acterized by the presence of aromatic clusters linked by 
aliphatic bridging groups. 

UV-induced dehydrogenation of the emerging 
arophatic cluster (or the parent a-C:H) introduces 
pentagonal rings, and hence curvature, in the evolving 
structure. 

The result is a large, vibrationally-excited fullerenic 
cage that shrinks down to the stable molecules Ceo and 
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C70 through emission of C2 groups. The very high en- 
ergy barrier for C2 ejection from Ceo prevents further 
dissociation. 

The young PNe where fullerenes have been observed 
appear as the optimised factories for the formation of 
these molecules. In their pre-PNe stages the formation 
of fullerene-like structure is not favoured because these 
objects lack sufficient UV photons to dehydrogenate and 
heat the dust. In their later stages the radiation field is 
so much harder that it will destroy H-rich carbonaceous 
particles before they can re-structure to form fullerenes. 
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